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Abstract 



The paper is devoted to prove the existence of a local solution of the Hamilton- 
CN ' Jacobi equation in field theory, whence the general solution of the field equa- 

tions can be obtained. The solution is adapted to the choice of the subman- 
ifold where the initial data of the field equations are assigned. Finally, a 
,S^ , technique to obtain the general solution of the field equations, starting from 

MH 1 the given initial manifold, is deduced. 
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Hamilton-Jacobi theory plays an important role within classical point particle 



I Introduction 
> 

O 

q^ mechanics, being a powerful tool to solve the evolution equations on the one hand, 

and the cornerstone for the transition from classical to nomrelativistic quantum 

f"*- ■ mechanics in Schroedinger' s formulation on the other. 

The formulation of Hamilton-Jacobi equation for field theories is quite old and 
possesses many different equivalent approaches [2]. It has been recently revised 
for its relationship with relativistic quantum mechanics [1] and for its geometric 
properties [12] . 

rS \ Like the classical case, the interest of the Hamilton-Jacobi method is placed 

in the possibility of using it to obtain solutions of the field equations. Indeed, the 
possibility of using the concept of complete integral to perform this operation is 
precluded, being its existence quite an exceptional case in field theory. 

In fact, while the equation itself is even less constraining than the mechanical 
one, being a single PDE in many unknown functions, an additional set of embed- 
dability conditions is needed to relate its solutions to the field equations. This 
results into a set of quite restricting conditions and the only available possibilities 
of finding solutions strongly depend on the nature of the particular problem under 
consideration. 



Many different authors pursued the aim of producing sets of particular solu- 
tions: some of them built single solutions starting from given solutions of the field 
equations [21 0], others use the concept of Baecklund transformation under particu- 
lar dimensional restrictions |5J. In a previous paper [1], a particular set of solutions 
was constructed; they were not constrained by dimensionality requirements, but 
the embeddability conditions resulted into a set of quite restricting compatibility 
relations on the initial data. 

This paper is devoted to proving the existence of solutions of the Hamilton- 
Jacobi equation in field theory, that allow to determine the general solution of 
the field equations. The only restrictions are placed in a set of quite general 
regularity conditions on the submanifold where the initial data are assigned. For 
this purpose, an algorithm is constructed, starting from a particular solution of 
the Hamilton-Jacobi equation, adapted to the choice of the initial manifold: this 
approach differs from the previous ones in its being sufficiently general to avoid 
that the embeddability conditions impose any restriction on the possible solutions 
of the field equations. 

Although the general solution was obtained, a complete integral, generating the 
whole parametric set of solutions, similar to the particle mechanics case cannot be 
constructed, being the solution still dependent on the choice of the initial surface. 
The solution of the field equation by means of the Hamilton-Jacobi method can 
however be obtained using an alternative concept of complete integral, adapted to 
the initial submanifold. This method is still a very partial algorithm, since it is 
not able to automatically include the embeddability conditions, that need to be 
integrated apart and is a current argument of research. 

The paper is organized as follows: 
Section II is devoted to recalling the main features of Hamiltonian field theory, 
which represents the background where Hamilton-Jacobi equation is placed. In 
particular the link between the solution of the field equations and the Hamilton- 
Jacobi theory is strictly revised. All the arguments are already contained in pQ. 
Section III describes the new technique for solving the Hamilton-Jacobi equation. 
The procedure builds a solution of the Hamilton-Jacobi equation, whence the 
solution of the field equations is deduced, for a general enough choice of the initial 
data. The proof of the existence of such a solution is limited to the class of analytic 
functions, being Cauchy-Kovalevskaya theorem needed to complete the argument. 
Section IV is devoted to showing possible future developments of the argument in 
the direction of defining a new notion of complete integral, adapted to the initial 
submanifold. 

In section V a simple example is given, in order to show how the whole machinery 
works. 



II Hamiltonian Field theory and the Hamilton-Jacobi 
equation 

The Hamilton-Jacobi theory represents a tecnique for solving the differential equa- 
tions of field theory in their Hamiltonian description. In this section, the most 
relevant aspects of Hamiltonian field theory will be briefly revised. The reader is 
referred to [1] for the details. The paper will directly deal with field equations in 
their Hamiltonian version; the relationship with the Lagrangian counterpart will 
be assumed as known. 

• Field theories are generally described from a geometric viewpoint taking a 
fiber bundle ir : E — > M into account, endowed with a set of local coordinates 
x 1 , . . . , x n , y 1 , . . . , y r , the first set denoting the basis coordinates on M and the 
remaining ones the local coordinates on the fibers. Within this context, a physical 
field is then represented as a section ip : M — > E, locally written in the form: 

y i = ip i {x^) i = l,...,r , /i = l,...,n. 

The base manifold M is supposed to be endowed with a (pseudo)-metric tensor $. 

• Hamiltonian field theory is developed within the framework provided by the 
bundle A n (E) of skew-simmetric n- forms on E (compare with p]). Let us take the 
bundle A™(E) of horizontal n- forms into account, whose elements are annihilated 
whenever one of its arguments is vertical. Then, let AV;(E) be the bundle of 1- 
contact forms on E, whose elements have the property of vanishing when at least 
two of their arguments are vertical vector fields. The elements of the two bundles 
can be locally described as follows: 

(j £ AJ(B) <^ a = po(a)ds , (la) 

a€A%(E) & a = p(a)ds + p^{cj)dy i AdSf, . (lb) 

Then, the bundles A™(E) and AV^(E) are respectively described by a system of 
local coordinates x tl ,y l ,pQ and x^,y l ,p,pf; moreover, A™(E) —>■ E is a vector 
sub-bundle of A%(E) -> E. The quotient bundle 11(E) := A^(E)/A^(E) will 
be henceforth called the phase space and endowed with a set of local coordinates 
x^,y l ,p^. It is easy to prove that A^iE) — > H(E) is an affine bundle, modelled on 
A?(£). 

• The assignment of the dynamical information is summarized giving a section 
S : p = -H{x^,y\p^) of the bundle A%(E) -» 11(E); the latter is called the 
Hamiltonian section and the function H is called the Hamiltonian function. It 
can be obtained from the knowledge of the Lagrangian of the theory through the 
Legendre transformation. The field equations, obtained from a variational problem 
on the phase space, can be written in terms of the Hamiltonian as follows: 

dy^_(M_ d^ = _dH_ 

dx^ dpf ' dx^ dy l 



The possibility of deducing eqs.([2|) from a Lagrangian formulation is related to the 
regularity condition 

, ( d 2 H \ , 

that will be tacitly assumed from now on. The calculation of the above determinant 
can be performed choosing an ordering, mapping the pair (/i, i) into a single index, 
running from 1 to n x r; the vanishing of the result is independent on the choice 
of the order. 

The field equations have the nature of a system of first-order partial differential 
equations for the functions y 1 = y l {x^),p^ = p^(x^). They are generally given 
together with a set of initial conditions, which essentially consist in the assign- 
ment of the fields and of their transverse derivatives over an (n — l)-dimensional 
submanifold i : N — >• M, locally expressed as 



^(z\ 



v n— 1\ 



\x = 1 . . . n 



where z , ■ ■ ■ , z n 1 represent a system of local coordinates on N. The whole set of 
initial data is therefore represented as: 



X' 



^{z A 
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is a vector field defined on N and transverse to it. Whenever 



the regularity condition ([3]) is satisfied, the equation 

.• dH 
dpj 
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^,y\pt) 



can be locally inverted with respect to pf. This allows to write eq.(j3]) in the form 
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-H{y{z A )^{z A ),${z A )) 



representing a (n — l)-dimensional submanifold of A^iE) 

• The Hamilton-Jacobi equation can be obtained determining a particular foliation 
of the Hamiltonian section S. For this purpose, let a £ A"~ (E) be a generic 
horizontal (n — l)-form, locally described as: 



a 



S a {x^y i )ds a 



(6) 



Its differential a G K^{Ej) is locally described by 

= — d^ — + — dt — y a ' ( } 

and represents a a section £ : E — > kQ(E), locally written as 

dS a (x^y*) a dS a (x»,y*) 

S: p= a* ' ^ = W ' () 

The latter is contained in the Hamiltonian section if the following condition is 
satisfied: 

^ +H U yi ^\^ . (9) 



dx a \ dy 1 

Eq. ([9|) is known as the Hamilton- Jacob i equation for field theories. 

• The solutions of the Hamilton-Jacobi equation are related to the solutions of 

field equations. Taking the derivative of ([9]) by y l we find that 

d dS a dH d 2 S a dH d ( dS a , „ ,, u ,.\ dH 

+ — - + 7: r- . — — = => - — ^--{x^,y l (x^)) = -T-- . 



cte a 9y l dy l dy % dyj dp^ dx a \ dy % J dy 

This means that the quantities pf = 4pr fulfill a subset of the field equations, 
for any choice of the functions y l (x fl ). Unluckily, the rest of the field equations 
are not necessarily fulfilled, and need to be checked by hand on the solution. The 
resulting condition is called embeddability condition and is represented as 

* BH '^,ym\ . (io) 



dx^ dpf \ ' ' dy 

A solution y % = y l {x v ), p^ = pf(x u ) of the field equations is said to be embeddable 
in a solution S^ = Sf l {x v ,y i ) of Q if condition CED holds. 

If fj, = 1, i.e. in the case of point particle mechanics, eq. f)10|) describes a velocity 
field on E for every solution of Q and its integrability is guaranteed by Cauchy 
theorem; therefore, it is possible to determine a family of solutions of the field 
equations embedded in every solution of the Hamilton-Jacobi equation, at least 
locally. In the general case, this is not necessarily true; the reader is referred to 
[1] for an analysis of the possible scenarios present in the literature. 

The approach adopted here will be that of determining a parametric set of solu- 
tions of the Hamilton-Jacobi equation, associated to a chosen initial submanifold 
N. This is made possible by the fact that eq. ^ has the nature of a single 
first-order partial differential equation in the n unknown functions S a (x^, y l ) and 
therefore possesses a great amount of equivalent solutions. This arbitrariness will 
be eventually fixed by the embeddability condition. The argument extends the 
analogous procedure presented in [I], where a particular choice of the candidate 



solution allowed to build particular classes of solutions, with the embeddability 
conditions giving rise to a set of compatibility requests on the initial data. This 
time, however, this choice will prove to be general enough to generate the whole 
set of solutions of the field equations, whenever the initial submanifold TV satisfies 
quite general regularity conditions. 

Ill The general solution 

Let us consider a field theory, where the Hamiltonian H is given, together with 
the surface N C M where the initial data are assigned. Then, let us consider a 
vector field X € D 1 (M), defined in an open subset U C M, containing the initial 
submanifold N, and locally written as 

X = X^(x\...,x n ) ° 



dx^ 



If X is chosen to be transverse to the initial submanifold N, its integral curves are 
locally determined as 

dr^ 

— = X»(xn -► x» = x»(Z,z A ) , (11) 

where £ is a parameter along the curves and z are the coordinates on JV, param- 
eterizing the starting point, corresponding to £ = 0. The set £, z , . . . , z n ~ l define 
a system of local coordinates around N, henceforth called adapted coordinate sys- 
tem, so that the submanifold N is represented by the slice £ = 0. Let us consider 
the candidate solution 

S^x^y i ) = ip(x",y i )X^ + Af(x^y i , (12) 

where Ai = A^{x^)-^jt is a family of r vector fields on M. The functions Af are 
arbitrary and the solution will depend parametrically on them. Substituting into 
(fT2|) we obtain 

For every choice of the functions Af, eq. (|13|) can be interpreted as a single partial 
differential equation for the function </?(£, y*). We notice that the coordinates 
z play the role of parameters, since no derivative is performed with respect to 
them. Moreover, the user-defined vector field X^ can be chosen in such a way that 
^jtt = 0. This choice allows to remove the dependence of the PDE on ip, thus 
simplifying the solution strategy. On the other hand, the choice is not coordinate 
independent and makes the following procedure strictly based on the chosen initial 
coordinates x , . . . x n . As a matter of fact, this limitation is not relevant and has 
to be thought on the same footing as the technique of separation of variables, that, 
when possible, is necessarily related to the choice of a suitable coordinate system. 
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Taking eq. (|ll|) into account, eq. (|13j) can be rewritten in the adapted coordinate 
system £,z as follows: 

Eq. (|14p is very similar to an ordinary Hamilton- Jacobi equation, where the role of 

the Hamiltonian is played by the function H + g-£- y l . 

In order to prove the existence of a solution and to determine the links with the 
solutions of the field equations, eq. (|14p will now be analyzed though the method 
of bicharacteristics (cmp. with |10j). For this purpose, given a single differential 
equation in a single unknown, we need to rewrite it as an algebraic expression of 
the form 

F(x i ,u,p i ) = , 

where x l are the independent variables, u represents the unknown function and 
Pi rs J^t. The solution is obtained integrating the ordinary differential equations 
describing the bicharacteristic curves, locally written as 

dx>__dF_ dp^^ dF_ dF_ du _ dF dF 

d£ dpi dS, dx l du % d£ du dpi 

and dragging the initial data along them. In the present case, taking the identifi- 
cations x l ~ (£, y l ), u ~ (p and pi ~ (no, Ui) into account, eq. (fl"4"|) can be rewritten 
as: 

u + ^y l + H(t,z A ,y\u l X» + A f ;) = . (16) 

Since eq. ()16p depends parametrically on the z s, the solution will be represented as 
a parametric family of functions. The differential equations for the bicharacteristics 
of eq.(fl"6|) take the form 



dtf_ dH 
d£ dp'. 



-^ = i^& z A , y\ UiX» + A£)X"(Z, z A ) (17a) 



^=- d S^ zA) -w { ^ AyUiXfl+Af) (17b) 

|f = -^z(Z, z A ,y\u t X» + A£)X»(Z, z A ) Ui (17c) 

and represents a system of 2n + 1 ordinary differential equations in the unknowns 

y*(Z,z A ),u(Z,z A ),u t (Z,z A ). 

First of all, we notice that eqs. lfTTh ). (fTTb ) are independent on u, so that eq. (fT7b ) 

can be integrated after solving the others. 

Moreover, the parametric dependence of ()17a ) and (|17b) on z makes their solution 

depend on 2n arbitrary integration functions of the parameters z , that can be 

eventually fixed assigning the initial data at £ = for the functions y l (0, z ) and 

u % {Q,z A ). 



The fact that the initial submanifold N is represented by the slice £ = 0, allows 
to set the integration functions from the knowledge of the initial data of the field 
equations ([!]). A first step is performed assigning y i (0, z ) = tp l (z ); the second is 
based on the application of the Legendre transformation, by means of which the 
normal derivatives of the fields at N are related to the values Uj(0, z ) by 



dy l 



dx* 1 



^_ iP . .1 ..... h .. ,„ .1. 



^V) = iTT^^V),^(0,Ar M + 4V • (18) 



By Dini's theorem, the possibility of calculating u,(0,z ) is guaranteed by the 
requirement 

M {im x * n ")*° ■ (19) 

The latter is not a direct consequence of the regularity condition on the Hamil- 
tonian function. Whenever the initial submanifold N is fixed, it should be thought 
as a restriction on the possible choices of the field X. 

The result of the above procedure consists in a candidate solution y % = y J (£, z ) 
fulfilling the initial data (JU) and still depending on the choice of the arbitrary func- 
tions Af. However, they are not solutions of the whole set of field equations, unless 
the embeddability conditions (fTOjl are satisfied. Taking eq. ffTTh ) into account,we 
notice that the subset of embeddability conditions 

*Lx* = ™x» 

is automatically fulfilled. Then, the embeddability conditions only consist in the 
system of equations 

f^f^V^ + Af)^ (20) 

that must be fulfilled by the functions y l (£,z ) and Uj(£, z ). The strategy of 
solution takes advantage of the complete arbitrariness of the functions Af and 
consists in imposing eqs. (|20p on them. 

For this purpose, let us decompose the fields Ai along the basis provided by 
the adapted coordinates as 

A.-A-— + k B —-(A~ + k Bdx ^ d 



'di * dz A ~\ l di l dz A J dx^ 
and let us write the embeddability condition in the form 

dy l dH ( ' A i * o <9z M \ dx^ 



^- 8i n^- 9 " ,( "- + ^ +fcf ^Ja? • <21> 



Whenever the regularity condition 



, / d 2 H dx 1 * dx u \ , . , 

det ( atfa^&W J * ° < 22 > 



is fulfilled, Dini's theorem allows to locally invert eq. (|2ip for the functions kf, 
namely 



.A _ , A I t A _i dy l 



kf = kf[Z,z A ,y\^ I ,u u A) . (23) 



where 



Substituting into (|17p we obtain: 

% = f(i,z A ,y\^UiX) (24a) 

-Q£=9i(€,z ,y\ui, t^, -fak, dz A dz B,Ai) (24b) 

- _^7* a * iV \ , dAj t dk£ (t A i &i_ dm ay 

ay 4 \ az^ / a£ az A az^ 1 az A oz A oz ti 

(25b) 
Eqs. (|24p represent a system of partial differential equations for the unknowns y l 
and Ui, in normal form relative to the derivatives with respect to £. 

We will now show that, for every choice of the function Ai, they possess a local 
solution around the initial manifold corresponding to £ = 0, represented by N. 
For this purpose we need to restrict the possible choices of the initial submanifold 
using the following regularity assumption. 

Definition 1 (Regular initial pair) A pair (N,X), where N is a submanifold 
of M having co-dimension 1 and X is a vector field of M transverse to N, is 
said to be a regular initial pair for the Hamiltonian H if satisfies the regularity 
conditions (|19j) and (|22p . together with 



det ( d2R X*X» - 92R Q2H X»X* 9XV 9XP A A A * (26) 

l [dp?d P » dp^dp" d P p r d P r dz A dz B ir J ^ u ' [b) 

where 

AB _ ( d 2 H dx^dx^ i' 1 
ij '~ ydpfdp" dz A dz B J 

We notice that the conditions defining the regular initial pairs are completely 
independent of the regularity condition on the Hamiltonian function but have to 
be considered as restrictions on the possible choices of the initial submanifold N 
and of the transverse vector field X. The hessian on the Hamiltonian plays the 
role of a metric tensor. 
We are now able to prove the following 



Theorem 2 Consider an analytic Hamiltonian section S : p + H(x^,y' t ,pf) = 
and a regular initial pair (N,X). Let B C S represent the initial data of the 
field equations (compare with ^) and suppose that X is analytic and satisfies the 
condition ^^ = 0. Then, the Hamilton- J acobi equation of H possesses a family 
of solutions of the form 

S^,y i ) = <p(x»,y i )X» + A?(x»)y i , (27) 

where Af are arbitrary functions on M . Moreover, there exists a choice of the 
functions Af such that the submanifold £ C A^iE), locally represented as 

S-l p^{i,z A ) = u t {i,z A )X^ + A>t{i,z A ) , 
.p=-H{^z A ,y\p») 

obtained dragging the initial the submanifold B along the bicharacteristic curves 
(|24p of the Hamilton- Jacobi equation, is a solution of the field equations for every 
given choice of the initial data. 

Proof We already showed that, substituting eq. fj27j) into (|14p and writing the 
equations for the bicharacteristics, we always determine a set of solutions of (|17|) 
matching all the possible initial data, and depending parametrically on Af. 
Being the regularity condition (|22p fulfilled, the parameters k A can be (locally) ob- 
tained from (|2ip by Dini's theorem, leaving us with a system of partial differential 
equations (|24p in normal form. Let us now define the functions 

F\^z A 1 y\u u y l A ,yl) :=yi- f {^z A ,y\y\,u i ) , (28) 



determined from (|24a ) choosing an arbitrary analytic function Ai. Applying Dini's 
theorem to (|28p it is possible to locally write the surface F l = in parametric form 



with respect to Ui whenever the condition det(^-) ^ is fulfilled. Taking eq. (|2"pp 
into account, it is easy to calculate 

df_\ = I d 2 H xuxv _ d 2 H d 2 H x »X x9xU 9xP A AB ) (29) 

duj J I dpfdp'j dp^dp v - dprdp^ dz A dz B ir J 

The non singularity of (|29p is guaranteed by eq. ()26p . We can therefore write 
eq. fl28|) as 

^ = m,Ay\^,^) ■ (30) 

We can also calculate the derivatives of (1301) as: 



di di + d V i 3^ + Qy* A dz A d£ + d y{ de [ ^ 
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(hn_ _ m dlhdyl dUj gy dUj 3 2 yi 

dz A dz A dyi dz A dy j B dz B dz A dy J d£dz A ' [ ' 

Substituting into eq. ([24"b ) we obtain: 



dUi dUi_(Hl_ dUi d2yj dU 1 d^_ f A { dy i d 2 y i d 2 y i 



<9£ <V d£ dy j A dz A d£ dyi d£ 2 ' *V dz A ' d£dz A ' dz B dz A 



being hi the right hand side of (|24b). after the substitutions of u« and its derivatives 
from (|30p and ()31 j) . This equation can be put into normal form with respect to 
whenever 



»V 



at 



2 



-§'"« 



^0 

This condition can be evaluated applying Dini's theorem to eq. ([2 
8U t __ fdF^' 1 dF k _ fdpY\k _ fdf^~ l 



dyi \du k J dyl \du k J 3 \du 



It is non singular by (|26p and ()29|) . This results into a system of second order 
partial differential equations in Cauchy-Kovalevskaya form as 

d 2 y i ( A t dy i d 2 y i d 2 y { 



de lr V ' ' dz A ' 0£0^ ' 0**0^ 

In the present case, this means that the system is in normal form for the second 
order derivatives with respect to £ and all other derivatives in any variable on the 
right hand side are of first or second order. Obviously, no second order derivative 
with respect to £ appears on the right hand side. Then, by a generalized version of 
Cauchy-Covalevskaya theorem |13j . whenever all the functions defining the system 
of partial differential equations are analytic, the latter possesses a local solution 
around the initial surface £ = 0, which represents the initial submanifold N. □ 

Now, taking eqs. ()17p into account, the general solution of the bicharacteristic 
equations (|24p can be written as 



y i = y l (tz A ,a i (z),(3 i (z)) , m = u i ^,z A ,a i (z),P i (z)) , (32) 

where the arbitrary functions a.i{z ) and f3 l (z ) can be calculated from the knowl- 
edge of the initial data of the field equations on N, using the same procedure 
reported in eq. (fTHj) . 



IV Further developments 

The results obtained up to now prove the existence of a solution of the Hamilton- 
Jacobi equation for every regular choice of the initial data of the field equations, 

11 



defined on a regular initial pair. From a computational viewpoint, the advantage 
of using the bicharacteristic equations over the usual field equation is not great. 
In this section we will show some ideas that can lead to a technique for solving 
the Hamilton-Jacobi equation using a direct approach. 

The argument uses the technique of the complete integral and is based on the 
fact that eq. (|14p is an ordinary Hamilton-Jacobi equation for the function <p(£, y z ), 
depending parametrically on the functions A%. For this reason, it possesses a 
complete integral of the form 

ip = <p(Z,z A ,y i ,a i ,A?) (33) 

satisfying the condition of essential dependence on the parameters 

The parameters ckj result to be first integrals of the problem, i.e. they are constant 
along the bicharacteristic curves. Taking eqs. (|17p into account, it is easy to prove 
the following 

Proposition 3 The functions j3 l = (3 l (£,z ) defined as 

F = ^(Z,z A ,y\a u A?) (35) 

don 

are constant along the bicharacteristic curves. 

Proof Taking the derivative along the bicharacteristic curves we obtain that 

dft _ d ( d<p\ dyi d 2 ^ 



d£ d£ \daij dti dyidcti 
Using Schwarz theorem and substituting from (|X4|) we obtain that 

_ a dA» dyi ay _ Q 

d£ doii dxt 1 d£ dyida 1 

The result comes from the fact that Af does not depend on on and 
d dH d 2 <p d 2 tp dy i 



da 1 dpf dy % da,i dy l dai d£ 

byeq. (fTTh). D 

As a consequence, the functions a, and /3 1 can be determined from the knowledge 
of the initial data. The essentiality condition (|34|) allows to local invert eq. (|35j) 
with respect to the y l , thus obtaining the candidate solutions of the field equations 
in the form: 

y i = y i ^,z A ,a i ,p i ,Af) . (36) 
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Eq. (|36p represents an actual solution of the field equations if it satisfies the em- 
beddability conditions. Taking the derivatives of (f55j) by z A , and substituting the 
embeddability conditions ()10p . we have that: 

dp _ d 2 y dH dx^ d 2 y daj d 2 ip 
dz A dctidyi dp^ dz A dctidctj dz A daidz A 

By Theorem 2, whenever the initial data take their values on a regular initial 
pair, there exists a choice of the function Af such that the functions (|36p are 
a solution of the field equations. However, solving eq.([37|) is as complicated as 
solving the field equations, making the whole argument still partial. Any further 
development of this research topic should therefore provide an algorithm including 
the embeddability condition within the solution technique of the Hamilton- Jacobi 
equation; this argument is still under analysis. 

V An example 

Let us consider the free scalar field in a 4-dimensional Minkowki spacetime M, 
endowed with a metric n = diag(l, 1, 1, —1). It Hamiltonian is given by 

H(x^y,p») = -rj^p u + -u 2 y 2 . 

Let us consider the problem of determining its evolution starting at N : x 4 = 
in the positive time direction. We choose X = -S-%, so that £ = x 4 . The surfaces 
£ = const are spanned by the coordinates z = x , A = 1..3) and are endowed 
with a positive definite metric 5ab- The regularity of the pair (N,X) can be 
proved observing that 

d 2 H _ 
dpvdp v ~ v ^ ' 

so that the conditions (fT9|) . ([22]) and ([26|) can be respectively written as: 



n^X^O , det( Vfll/ S^S u B )^0 , n lu ,Xi t X v £Q . 

This shows that a solution can be obtained whenever X is a timelike vector field 
and the pullback of the metric tensor on the initial surface is non singular. 
The Hamilton-Jacobi equation becomes 

dS» 1 dS^dS" 1 o 2 

d^ + 2^^y-^y7 + r y = ° ■ 

Substituting (|27|) we have that 



8i + fe?» + 2"- M " [fy) + 2^ A " A " + ^ X 8y A " + 5" * 
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We notice that r] )lu X ,l X 1 ' = — 1 and, choosing A = A 4 = 0, we also have that 
•q iaf X' i A v = 0. Substituting h A := A A ,A = 1..3 and \\h\\ 2 := i]^A^A u = 
$ABh h > 0, we have that 



2 



dip dh A 1 (dp\ 1 2 , 1 2 2 



Writing ui ~ ■# , the equations of the bicharacteristics are 
dy du\ dh A 

d~r~ ui ; ^ = ~ fiy ~d^ • 

Its separation easily leads to 

dy d 2 y AB _ 2 2 



9£ dz A dz B 



5 A " = mV 



that are identical to the field equations for the free scalar field. 

In order to show what difficulties rise in the application of the arguments of 
Section IV, a complete integral will be determined using the following strategy. 
Let us consider a polynomial solution 

V = -tiy 2 +a(£,z A )y + b(Z,z A ) . 

Substituting into (|3"5j) we have that 

da db 1 n 1 I|7||2 dh A 

ez* + Tr2 a aw+ 2 m + d^ y = ■ 

Separating the linear and the homogeneous term in y we can solve for a(£, z ) and 
b(£,z A ) as 

/dh A 



b^z A ) = \ j a 2 <% + \ j \\h\\ 2 <% 



where the integration constant a(z ) is the required parameter for the complete 
integral. We notice that the dependence on a is essential, since 



*v __ ^ 



dyda 
Then we can calculate 

dip 



e« 7^0 



9a y 
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whence 

y(f , z) = e-^P - e~^ I a e^d£ . (39) 



Eq. pUj) represents the candidate solution of the problem. It also has to fulfill the 
embeddability conditions, that can be calculated as follows: 

d(3 _ u , , e 2 ^ da , _ 2uf f __ u , d 2 h 



dz A 



^+wW^w« ' 



The difficulties of this last equation , representing the embeddability conditions, 
are similar to those of solving the field equations directly. In fact, it can be easily 
written in the form 

d 2 h B _ d 2 h A _ 2 . 
dz B dz A d£> 2 -V n * , 

whence the solution can be easily found. 
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